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Slowly Pulsating B and /3 Cephei are k, mechanism driven pulsating B stars. That n mechanism works since a peak in 
the opacity due to a high number of atomic transitions from iron-group elements occurs in the area of logT « 5.3. 
Theoretical results predict very few SPBs and no /3 Cep to be encountered in low metallicity environments such as the 
Small Magellanic Cloud. However recent variability surveys of B stars in the SMC reported the detection of a significant 
number of SPB and j3 Cep candidates. Though the iron content plays a major role in the excitation of /3 Cep and SPB 
pulsations, the chemical mixture representative of the SMC B stars such as recently derived does not leave room for a 
significant increase of the iron abundance in these stars. Whilst abundance of iron-group elements seems reliable, is the 
opacity in the iron-group elements bump underestimated? We determine how the opacity profile in B-type stars should 
change to excite SPB and jS Cep pulsations in early-type stars of the SMC. 
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1 Introduction 

The Small Magellanic Cloud is a dwarf galaxy satellite of 
our galaxy, that presents very low metallicity: Z w 0.001 — 
0.004 (Buchler 2008). Recent theoretical study (Miglio et 
al. 2007a) employing the updated solar chemical mixture 
AGS05 (Asplund, Grevesse & Sauval 2005) find that at Z= 
0.005, no excitation of /3 Cep modes and a few excited SPB 
modes in B stars are expected. Therefore SPBs are expected 
to be rare while (3 Cep are thought to be absent in the SMC. 
However, recent detections of B-type pulsator candidates in 
the SMC (e.g. Karoff et al. 2008; Diago et al. 2008; Sarro 
et al. 2009) are now challenging our understanding of the 
main-sequence B-type pulsators. 

In a previous work (Salmon et al. 2009), we tried to ex- 
plain the existence of these B-type pulsators in the SMC by 
considering a chemical mixture representative of the SMC 
B stars. From a synthesis of the literature, we derived it and 
employed it to compute new theoretical models. Our results 
indicated that the minimal values of Z to obtain excitation 
respectively of SPB and /3 Cep are both over the mean Z 
value of SMC B stars, even when the maximal error box 
value is taken into account. 

We foresee another possibility to reconcile observation 
and theory. The mechanism of excitation in /? Cep remained 
for long unknown: Simon (1982) suggested that a revision 
of opacity values for heavy elements by a factor 2 or 3 could 
explain excitation in /? Cep stars that was unexplained at that 
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time. Several years later, new computations of opacity by 
OPAL group (Rogers & Iglesias 1992) led to an increase by 
a factor 3 of the opacity due to iron-group elements, hence 
resolving B-type pulsators problem. More recently, an un- 
derestimation of current stellar opacity values has already 
been invoked to solve other remaining problems such as the 
excitation of low and high frequencies in galactic /3 Cep 
pulsators, the discrepancy between the solar model with the 
updated solar mixture and helioseismic data and the model- 
ing of the slowly-rotating (5-Scuti 44 tau (for a review, see 
Montalban & Miglio 2008). 

These studies raise the following question: is stellar opac-| 
ity still underestimated nowadays? In this work we estimate 
by which factor the opacity should be increased to excite 
pulsations in SMC B stars. We proceed by modifying the 
opacity values used in the theoretical computations. 

We briefly recall our previous results in Section 1. In 
Section 2, we explain the way we proceed to modify the 
opacity in our models. We then present the results of our 
non-adiabatical computations in Section 3 and the conse- 
quence of an increase of the opacity on the excitation of 
B-type pulsators. 

2 Effect of the chemical mixture 

Motivated by the apparent inconsistency between observa- 
tion and theory for B-type pulsators in the SMC, we ana- 
lyzed the effect of considering a chemical mixture represen- 
tative of SMC B stars (Salmon et al. 2009). After a careful 
and critical review of the literature of SMC age-metallicity 



®WILEY 

3 iTiterScience' 



© 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Wcinheim 



790 



T.H.E. Editor & G.H. Ostwriter: Instructions for authors 




Z=0.0024 M=12M. SMC mixture 





-1 1 1 1 1 1 1 1 1 1 1 1 


' 1 ' ' - 




Z=0.0024 <cxl 






Z-0.003 «Kl 


'■■. 


0.5 


Z=CI.004 KKl / , 

Z=0.005 KXl / / 


"\^; 




Z=0-01 KxX I jlu 


•■■.:::r' 




Z=0.Q024 /cxl.30 / |.7 


■ 





/ 


- 


0.5 




. 1 . . 



6 

loE T 



log T 



Fig. 1 Left panel : profiles of the enhanced opacities for a given structure {p and T given for a YIMq, Z=0.0024,Y=0.28 
on ZAMS model). The greater is the A factor, the wider is the opacity peak. Right panel : same as to the left, except that 
the 30% increased opacity at Zs„c is compared to opacities at different Z, but not parametrically modified. 



relation on the one hand and of chemical abundances deter- 
mination of young objects (B stars and cool supergiants) on 
the other hand, we determined two main features for young 
SMC objects: 

- one expects [Fe/H]=-0.5±0.2 with respect to AGS05 so- 
lar chemical mixture for the youngest stars of the SMC 
(Carrera et al. 2008 and ref. therein) due to their age and 
localization: the SMC B stars benefit from the largest 
chemical enrichment among all the stars of the galaxy. 

- the chemical mixture representative of the SMC B stars 
(see Salmon et al. 2009), implying Zsmc=0.0024±0.0006 
with [Fe/H]=-0.7±0.2. The iron content we derived is 
in agreement with the one based on age-metallicity re- 
lations previously quoted. 

We thus computed models with the stellar evolutionary 
code CLES (Code Liegeois d'Evolution Stellaire; Scufflaire 
et al. 2008), using OP opacities determined for the SMC 
B stars chemical mixture. The non-adiabatical analysis was 
performed with the MAD code (Dupret 2001). 

For models with masses ranging from 2.5 to 16 M©, we 
did not obtain excited modes for Zsmc- Moreover, we did the 
same analysis on models of 4 and 12 M©, respectively rep- 
resentative of SPB and /? Cep stars, for different values of 
Z. The results led to the conclusion that Z should be greater 
than 0.004 (corresponding to [Fe/H]^-0.45) to obtain ex- 
cited SPB-type modes and greater than or equal to 0.007 
([Fe/H]~-0.2) to obtain excited /3 Cep-type modes. 

3 Parametric increase of the opacity 

As we cannot explain the existence of /3 Cep pulsators (SPB 
being marginally compatible) in the SMC employing Zsmc 
and a representative chemical mixture in our models, the 
following question follows: is the stellar opacity underes- 
timated? Recent studies explored the effect of increasing 



opacity by enhancing the iron abundance in stellar models 
(see Ausseloos et al. 2004; Pamyatnykh, Handler & Dziem- 
bowski 2004; Miglio et al. 2007b). In the same idea we pro- 
ceed to an ad hoc and local change in the opacity used in 
our models, in order to estimate which enhancement of the 
opacity could lead to excitation of SPB- and /3 Cep-type 
modes at Zsmc- Our way of proceeding allows us to empha- 
size directly the consequences of a modification in the em- 
ployed fundamental physics. 

Iron is the element mainly contributing (sa 67%) to the 
peak of opacity in the layers with T w 200, 000 K, responsi- 
ble for excitation in B stars. Thus we decided to consider an 
increase of the opacity which would be due to an underes- 
timation of the iron contribution. We proceed by modifying 
the opacity (OP one, calculated for the SMC B stars mix- 
ture) as follows: 



K'(r) =k(T) X <^ 1 + ^ exp 






(1) 



where k is the Rosseland mean opacity and Tp^ and ope 
are respectively equal to 5.30 and 0.5. These values have 
been chosen to match accurately the contribution of iron to 
the opacity in the area of T sa 200, OOOi^. The factor of 
increase. A, has been set to different values, namely 0.05, 
0.10, 0.20, 0.30, 0.50 and 1. The modified opacities are 
represented in Fig. [T] (left panel). The modified iron peak 
opacity value is equivalent to a case between Z=0.004 and 
^=0.005 for unchanged opacities, as illustrated in Fig. 1 
(right panel). However the modified profile has a different 
slope in the iron bump area than in the case of unchanged 
profiles, this slope impacting on the damping of modes. The 
excitation depends not only on the value of opacity near 
log T=5.3 but also on its change with T. 
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Fig. 2 The 77 parameter of modes (reported to the logarithm of their period) for a typical SPB (resp. /3 Cep) model of the 
SMC in the left (resp. right) panel. The 77 parameter indicates the rate of growth of a mode. A positive value (resp. negative) 
indicates that the mode is excited (resp. damped). 



Table 1 Value of the increasing factor A necessary to 
obtain excitation of p or g modes is reported for the different 
combinations of model parameters (Z and M). 
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4 Stability analysis for models with modified 
opacity 

We calculate models of 4 and 12 Mq, using the opacities 
which we modified for the different factors of increase A 
and Z=0.0024, 0.003, 0.004, 0.005 and 0.01. In Fig. [3] we 
present in a HR diagram the main sequence evolutionary 
tracks of models for masses between 2.5 and 16 Mq, adopt- 
ing unchanged opacity. On the same figure are represented 
the 4 and 12 Mq models at Zshc computed with the dif- 
ferent increased opacities. For both masses, the position of 
main sequence tracks in the HR diagram is not considerably 
affected. 

Results of the non-adiabatic calculations are shown for 
the 4 and 12 Mq models at Zsmc in the left and right panels 
of Fig.l2]respectively. It appears that SPB modes are excited 
for an opacity multiplied by 2 {A=l), while no /? Cep modes 
are excited for the factors of increase A we used. In Fig. 4 
(left panel), one remarks that excited g modes are found at 
the end of the evolution on the main sequence for the SPB 
model in which the iron peak is increased by a factor 2. 
Excited frequencies of the 4:Mq model at Zsmc and for an 
increase A=l are compared to the observational frequencies 
of the detected SMC SPBs in Fig. 4 (right panel). 




log Tj, 



Fig. 3 Evolutionary tracks for non-modified and modi- 
fied opacity models in the HR diagram. Indicative instability 
strips of SPB and (3 Cep-type pulsations are obtained from 
models computed with Z=Q.Ql and k x 1. 



Finally, we report in Table [T] the minimal factor of in- 
crease A that excites p or g mode in the 4 or 12 Mq models 
for the metallicities we considered. We find excited g modes 
in the case Z=0.0024, M=4 M© and A=l. When consider- 
ing Z=0.003 (upper limit within the l-a error box of the 
Z of the SMC B stars), g modes are excited in the 4 M0 
model for A=0.50. A slight increase (10%) of the opacity 
at Z=0.004 excites g modes in the 4 Mq model. Excited p 
modes are only found in the /3 Cep model (12 M0), appear- 
ing at Z=0.005 for a factor 2 {A=l) of increase in the opac- 
ity. For the same conditions, excited high-order g modes are 
also obtained. 
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Fig. 4 Left panel : Evolutionary tracks for a 4 M© model computed for several increases in the opacity Iron peak. Models 
with theoretically excited p or g modes are marked by a pentagon or a triangle. Right panel : on the left part of the figure 
is represented the logarithm of the frequency (in c/d) for SPB modes observed in the SMC (Diago et al. 08; Karoff et al. 
08 and Sarro et al. ). On the right part is reported the logarithm of the frequency to the log T^/f of the model in which the 
mode is found excited. 



Under the gaussian profile assumption that we have made| 
to modify the opacity profile, one remarks that increasing 
the opacity in the iron peak does not change the intrinsic 
nature of modes excited in our models: SPB and /3 Cep rep- 
resentative models remain of the same nature. In the 12 Mq 
model at Z=0.005 and A=l as excited low order p modes 
appear, so do high order g modes. From an observational 
point of view, considering opacities as underestimated by 
more than 50%, SPB modes should be detected in SMC B 
stars with Zs„c=0.0024 ([Fe/H]=-0.7) while /3 Cep modes 
should be observed in the stars with Z=0.005 ([Fe/H]=- 
0.35), i.e. a metallicity within 5-a error box on Zs„c- 



and its consequence on the pulsation driving efficiency will 
be refined in the future. 

The stellar depth where stands the opacity peak also 
plays a role both on the damping of modes and on the ex- 
cited frequencies spectrum in B stars. Other elements of the 
iron-group (Ni, Cr and Mn), contribute to that bump in the 
opacity at slightly different temperatures affecting position 
of the opacity peak. If we suppose an underestimation of 
the iron opacity, it would be reasonable to consider the same 
for the other iron-group elements opacity. Thus we will take 
them into account in the next step of our theoretical model- 
ing of SMC B stars. 



5 Conclusion 

We have investigated the effect on the excitation of modes 
in SMC B stars of an underestimation of stellar opacity val- 
ues in the peak due to iron. We proceed by multiplying the 
opacity used in our model computations by a gaussian pro- 
file centered on the temperature where contribution of iron 
to opacity is maximum, i.e. logT=5.3. The results indicate 
that at the metalUcity of the SMC B stars, Zs„c=0.0024, SPB 
modes could be excited for an augmentation greater than 
50% of the opacity. No /3 Cep modes are excited at Zjmc for 
an increase from 5% to 100% of the opacity. Nevertheless, 
(3 Cep modes are found excited for an increase greater than 
50% but at Z=0.005, that is 5-az (a;j=0.0006) away from 
the metallicity we derived for SMC B stars. These values 
of increase must be taken with care since they are obtained 
for opacity modified according to an ad hoc defined math- 
ematical function that determines the slope of the opacity 
profile. The impact of the choice of mathematical function 
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